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A scheme for entangling atoms (ions) located close to or encapsulated inside a carbon nanotube 
is investigated using the photon Green function formalism for quantizing electromagnetic fields in 
the presence of quasi-one-dimensional absorbing and dispersing media. Small-diameter metallic 
nanotubes are shown to result in a high degree of the quantum bit entanglement for sufficiently long 
times. 
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In spite of impressive experimental demonstrations of 
basic quantum information effects in a number of differ- 
ent mesoscopic systems, such as quantum dots in semi- 
conductor microcavities, cold ions in traps, nuclear spin 
systems, atoms in optical resonators, Josephson junc- 
tions, etc., their concrete implementation is still at the 
proof-of-principle stage (see, e.g., Ref. Q for a review). 
The development of materials that may host quantum 
coherent states with long coherence lifetimes is a criti- 
cal research problem for the nearest future Q ■ There is a 
need for the fabrication of quantum bits (qubits) with co- 
herence lifetimes at least three-four orders of magnitude 
longer than it takes to perform a bit flip. This would 
involve entangling operations, followed by the nearest 
neighbor interaction over short distances and quantum 
information transfer over longer distances. It is thus of 
vital importance to pursue a variety of different strategies 
and approaches towards physically implementing novel 
non-trivial applications in modern nanotechnology. 

In this Letter, we study a scheme for entangling atoms 
(ions) located close to or encapsulated inside a car- 
bon nanotube (CN). CNs are graphene sheets rolled- 
up into cylinders of approximately one nanometer in di- 
ameter. Extensive work carried out worldwide in recent 
years has revealed intriguing physical properties of these 
novel molecular scale wires Nanotubes have been 
shown to be useful for miniaturized electronic, mechani- 
cal, electromechanical, chemical and scanning probe de- 
vices and as materials for macroscopic composites 0, . 
Recent experiments on the encapsulation of single atoms 
(ions) into single-walled CNs i[5j and their intercalation 
into single- wall CN bundles [6j, along with the progress 
in growth techniques of centimeter-long small-diameter 
single- walled CNs Q, stimulate the study of dynamic 
quantum coherent processes in atomically doped CNs. 

It was shown recently that the relative density of pho- 
tonic states (DOS) and the atom-vacuum-field coupling, 
respectively, near CNs effectively increase due to the 
presence of additional surface photonic modes coupled 
to CN electronic quasiparticle excitations ||. In small- 
diameter CNs the strong atom-field coupling may oc- 
cur 0, 0, — the property that is known to facili- 
tate the entanglement of spatially separated qubits [l^ . 



13] . Qualitatively, in terms of the cavity quantum elec- 
trodynamics (QED), the coupling constant of an atom 
(modelled by a two-level system with the transition 
dipole moment (La and frequency uja) to a vacuum field 
is given by fig = {2nd A Tiuj A /V) 1 / 2 with V being the 
effective volume of the field mode the atom interacts 
with (see, e.g., Ref. 0]). For the atom (ion) encap- 
sulated into the CN of radius R cn , V ~ irR 2 n (\A/2) 
that is ~ 10 2 nm 3 for CNs with diameters ~ 1 nm 
in the optical range of Xa ~ 600 nm. Approximating 
d,A~er~e(e 2 /htjJA) [15|. one obtains %~0.3 eV. On the 
other hand, the "cavity" linewidth is given for uja in res- 
onance with the cavity mode by h~f c = 6irhc 3 / 'uj^^uj a)V , 
where £ is the atomic spontaneous emission enhance- 
ment /dehancement (Purcell) factor ^ij. Taking into ac- 
count large Purcell factors ~ 10 7 close to CNs Q, one 
arrives at ft/Y c ~0.03 eV for 1 nm-diameter CNs in the op- 
tical spectral range. Thus, for the atoms (ions) encapsu- 
lated into small-diameter CNs the strong atom-field cou- 
pling condition g/"f c 2> 1 is supposed to be satisfied, giv- 
ing rise to the rearrangement (" dressing" ) of atomic lev- 
els and formation of atomic quasi-one-dimensional (ID) 
cavity polaritons 0, 0, 0. The latter ones are simi- 
lar to quasi-OD excitonic polaritons in quantum dots in 
semiconductor microcavities 16], that are currently be- 
ing considered a possible way to produce the excitonic 
qubit entanglement [l7j . Here we suggest an alternative 
way to generate the qubit entanglement by using quasi- 
1D atomic polariton states in CNs. We show that small- 
diameter metallic nanotubes result in sizable amounts 
of the two-qubit atomic entanglement. Envisaged appli- 
cations of our scheme range from quantum information 
transfer over long distances (centimeter- long distances, 
as a matter of fact, since centimeter-long small-diameter 
single- walled CNs are now technologically available ) to 
novel sources of coherent light emitted by dopant atoms 
in CNs. 

We use our previously developed photon Green func- 
tion formalism for quantizing electromagnetic field close 
to quasi- ID absorbing and dispersing media 0,0] . Rep- 
resenting such a medium, the (achiral) CN is considered 
to be an infinitely long, infinitely thin, anisotropically 
conducting cylinder. Its (axial) surface conductivity is 
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taken to be that given by the 7r-electron band structure in 
the tight-binding approximation with the azimuthal elec- 
tron momentum quantization and axial electron momen- 
tum relaxation taken into account. Two identical two- 
level atoms, A and B, are supposed to be positioned at 
their respective equivalent places ta,b close to the CN. 
We assign the orthonormal cylindric basis {e r , e tpi e z } in 
such a way that e z is directed along the nanotube axis 
and, without loss of generality, ta — ta^t = {^,0,0}, 
vb = {fBt Wb, zb}- The atoms interact with a quantum 
vacuum electromagnetic field via their transition dipole 
moments that are assumed to be directed along the CN 
axis, &a,b = d z e z . The contribution of the transverse 
dipole moment orientations is suppressed because of the 
strong depolarization of the transverse field in an iso- 
lated CN (the so-called dipole antenna effect 18] ). We 
also assume the atoms to be located far enough from 
each other, to simplify the problem by ignoring the in- 
teratomic Coulomb interaction. The total secondl y qu an- 
tized Hamiltonian of the system is then given by 



H = 



TiCbi 



duo huo / dR /t(R, w)/(R, w) + J2 ~« 1 a 



i=A,B 



/'L/dR[gW(r i ,R,w)*t 

i=A,B J ° J 

g(-)(r 4 ,R,w) a- i ]/(R,w)+h.c, 
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where the three items represent the medium-assisted 
(modified by the presence of the CN) electromagnetic 
field, the two-level atoms, and their interaction with 
the medium-assisted field, respectively. The operators 
f'(R,oj) and /(R, oj) are the scalar bosonic field oper- 
ators defined on the CN surface assigned by the radius- 
vector R ={R cn , <f>, Z} with R cn being the radius of the 
CN. These operators create and annihilate the single- 
quantum bosonic-type electromagnetic medium excita- 
tion of the frequency u> at the point R of the CN sur- 
face. The Pauli operators, &i = \L)(ui\, u\ = \ui)(L\ 
and ai z — \ui){ui\ — \L)(L\ with i = A or B, describe 
the atomic subsystem and electric dipole transitions be- 
tween its two states, upper \m) with either of the two 
atoms in its upper state and lower \L) with both atoms 
in their lower states, separated by the transition fre- 
quency to a- This (bare) frequency is modified by the 
diamagnetic (~A 2 ) atom- field interaction yielding the 
new renormalized transition frequency uoi = uoa[1- — 
2/(%uoa) 2 Jo rfw/dR|g ± (ri, R, uo)\ 2 ] in the second term of 
the Hamiltonian. The dipole atom-field interaction ma- 
trix elements are given by = g^ ± (lu/uja)^ where 
g^") = -i{Auo A /c 2 )^i:Tiuo Re a zz (uo) d/ (ll) G zz (r 4 , R, w) 
with -^^Gzz^i, R, uo) being the zz-component of the 
transverse (longitudinal) Green tensor (with respect to 
the first variable) of the electromagnetic subsystem and 
<7 zz (w) representing the CN surface axial conductiv- 
ity. The matrix elements g-^H) have the property of 



/dR|g x (ID(r 2 ,R,u;)| 2 - {^^{uoA/oofT^uo)^^^) 
with ^\r 2 ,uj) = lm ± ^Gzh (]l \r l ,r l ,Lu)/lmG°zz(uj) be- 
ing the transverse (longitudinal) local photonic DOS 
functions and T (uj) = 8Truo 2 d 2 ImG zz (w)/3?ic 2 represent- 
ing the atomic spontaneous decay rate in vacuum where 
\m.G zz {oS) ~lu/6ttc is the vacuum imaginary Green ten- 
sor zz-componcnt. The Hamiltonian Q involves only two 
standard approximations: the electric dipole approxima- 
tion and the two-level approximation. The rotating wave 
approximation commonly used is not applied, and the 
diamagnetic term of the atom-field interaction is not ne- 
glected (as opposed to, e.g., Refs. H,E|). 

For single-quantum excitations, the time-dependent 
wave function of the whole system can be written as 



4(Qi— oj)t I 



«>|{o» 
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i=A.B 



+ J duo jdR C L (R, w, t) e -*b>-»)* \ L) | 1 (R, uo)) , 

where uo = 5Zj=AB^/^' K^}) ls the vacuum state of 
the field subsystem, |{1(R, oj)}) is its excited state with 
the field being in the single-quantum Fock state, Cm and 
Cl are the respective probability amplitudes of the upper 
states and the lower state of the system. For the following 
it is convenient to introduce the new variables C±(t) = 
[Cua (t) i C UB (t)]/ \[2 that are the expansion coefficients 
of the wave function (J2J in terms of the maximally en- 
tangled 2-qubit atomic states |±) = (\ua) ± \ub))/V2- In 
view of Eqs. (Q) and 101, the time-dependent Schrodinger 
equation yields then 



C±(r) 



dr'K ± (T-T')C ± (T') + f ± (r), (3) 



where 



K±(t -t') = -I dx^^irA, v B ,x)e-^-^ T - T '\ 



2tt 



(4) 



e(r A ,r B ,x) = -± [Z x (r A ,x)±Z x (r A ,r B ,x)] (5) 
+ (H r A,x) ±^(r A ,r B ,x), 
Mr) =--j= fdr'K ± (r - r') C Ua (t'), (6) 

and we have, for convenience, introduced the dimension- 
less variables To = ^0/270, x = huo/2^, and r = 2"fot/h 
with 7o = 2.7 eV being the carbon nearest neighbor hop- 
ping integral (ft/270 = 1.22 xl0~ 16 s) appearing in the CN 
surface axial conductivity a zz . The functions /±(t) are 
only unequal to zero when the two atoms are initially in 
their ground states, with the initial excitation residing in 
the nanotube. Eq. © assumes that this situation is real- 
ized by selecting the time origin to be right after the time 
interval At, that is necessary for the (excited) atom A to 
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decay completely into the nanotube photonic modes, has 
elapsed. The two-particle local photonic DOS functions 
t x W{T A ,T B ,x) =Im M ^G^ l) (r A ,r B ,x)/lmG° z (x) are 
the generalizations of the DOS functions (r A , x) (see 

Refs. [H EH). They are of the form of £ x = 1 + 

C" = C" with the medium-dependent contributions given 

for r A , B < R cn by 

U ('A,r fl ,aj)J t 7 C (J/ « J 

(7) 

if 2 [vu (R cn )x]I p [vu(r A )x]Ip[vu(r B )x] cos [u (z B ) xy] 
1 + s(R cni x)v 2 I p [vu(Rcn)x]Kp[vu(R cn )x] 

where I p and K p are the modified cylindric Bessel func- 
tions, v(y) = v/y 2 - 1 , u(r) = 2jor/hc, and s(R cn ,x) = 
2iau(R cn )xa zz (R cn ,x) with cf zz = 2irh<r zz /e 2 being the 
dimensionless CN surface conductivity per unit length 
and a = e 2 /Tic = 1/137 representing the fine-structure 
constant. The integration contour C goes along the real 
axis of the complex plane and envelopes the branch points 
y = ±l of the function v(y) in the integrands from below 
and from above, respectively. For r A ,B > Rcn, Eq. Q is 
modified by the replacement r A ,B <~* Rcn in the Bessel 
function arguments in the numerator of the integrand. 

The entanglement of two quantum bits occurs when 
the 2-qubit wave function cannot be represented as 
a product of the two 1-qubit states in any basis. To de- 
termine this quantity in our particular case, we follow 
the recipe based on the "spin flip" transformation and 
valid for an arbitrary number of qubits (see Ref. flif 
for details). First, we define the reduced density ma- 
trix p A b(t) = \iPab(t)){i/iab(t)\ = Trf ie i d \ip(r))(ip(T)\ 
describing the bipartite atomic subsystem in terms of the 
wave function J2J of the whole system. Next, we intro- 
duce the "concurrence" G(^ab) = ab\iI> ab)\ where 

IVU-b) = &y°'y\' t P A B) wu; h o'y^ being the Pauli ma- 
trix that represents the " spin flip" transformation in the 
atom A(B) single-qubit space. This, after some alge- 
bra, becomes C[^ AB {r)] = |C 2 (r) - C 2 (r)| with C±(r) 
given by the integral equation (J3J. Finally, the degree of 
the entanglement of the 2-qubit atomic state \i])ab) is 
given by E[ip AB (r)] = h{l + y/1 - C[^ab(t)] 2 /2} where 
h(y) = -ylog 2 y- (1 - y) log 2 (l - y). 

The entanglement E[^as(t)] is maximal when the co- 
efficients C+ (r) and C_ (r) are maximally different from 
each other. For this to occur, the functions £^(ta,Tb,x) 
in Eqs. should be different in their values. These 

are determined by £ ^\r A ,r B ,x) whose frequency be- 
havior is determined by the CN surface axial conduc- 
tivity a zz . We computed £ from Eq.© with a zz 
calculated beforehand in the relaxation-time approxima- 
tion (relaxation time 3 x 10 -12 s at temperature 
300 K HHQ. Then, the integral equation was solved 




FIG. 1: (Color online) (a) Transverse local photonic DOS 
functions £ (a;) for the two- level atom in the centers of the 
four 'zigzag' nanotubes. (b) Normalized two-particle local 
photonic DOS functions £ + (solid lines) and £ - (dashed lines) 
taken at the peak frequencies of £ x (a;) [see (a)], as functions 
of the distances between the two atoms on the axes of the 
(10,0) (lines 1; x = 0.29), (11,0) (lines 2; x = 0.25) and (12,0) 
(lines 3; x = 0.24) CNs. (c) Same as in (b) for the two atoms 
on the axis of the (9,0) CN (a; = 0.32). (d) Two-particle DOS 
functions £ + (x) (solid lines) and £,~(x) (dotted lines) for the 
two atoms located in the center of the (9,0) CN and separated 
from each other by the distances of 2.1R cn « 7.4 A (lines 1), 
4.2i? cn «14.8A (lines 2) and 6.3i? c „ ^22.2 A (lines 3) [shown 
by the vertical lines in (c)]; the dashed line shows £ x (a;) for 
the atom in the center of the (9,0) CN. 



numerically to obtain C±(t) and E^ab ("?")]• The vac- 
uum spontaneous decay rate was estimated from the 
expression Tq(x) ~ o?x valid for atomic systems with 
Coulomb interaction [l5j |. The atomic transition fre- 
quency x A was assumed to be the local DOS reso- 
nance frequency. The simplest, most interesting case was 
considered where both of the atoms are positioned in the 
center of the CN. 

Shown in Fig. ^ (a) is a typical frequency behavior of 
the one-particle transverse photonic DOS ^(r^ = 0,a;) 
in the infrared and visible frequency range x < 0.4 for 
the atom inside 'zigzag' CNs of increasing radii. The 
DOS resonances are seen to be much sharper for metal- 
lic CNs (m = 3q, q = 1,2, ...) than for semiconducting 
(to =/= 3q) in agreement with the fact that this frequency 
range is dominated by the classical Drude-type conduc- 
tivity which is larger in metallic CNs ||| . Figure ^ (b) 
shows the normalized two-particle local photonic DOS 
functions £ x taken at the resonance frequencies of the 
respective ^'s [Fig. 1 (a)], as functions of the distance 
between the two atoms on the axes of the (10,0), (11,0) 
and (12,0) CNs. The values of £ + and £~ are seen to be 
substantially different from each other before they reach 
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FIG. 2: (Color online) (a),(b),(c) Upper decay of the ini- 
tially excited atom A (lines 1) and initially unexcited atom B 
(lines 2), and the entanglement (lines 3), as functions of di- 
mensionless time for the two atoms located in the center of 
the (9,0) CN and separated from each other by the distances 
of 2.1R cn » 7.4 A, 4.2R cn « 14.8 A and 6.3R cn » 22.2 A, re- 
spectively [the situation shown in terms of local photonic 
DOS in Figs. Uc),(d)]. (d) Short-time and long-time entan- 
glement evolution of the initially fully entangled atoms sepa- 
rated by 7.4 A; lines 1 and 2 are for C+(0) = 1, C_(0) = and 
C+ (0) = 0, C- (0) = 1, respectively. 

their limit values of £+ = £~ = £ x + £H » 2£ x at suf- 
ficiently large interatomic separations. For metallic CNs 
the functions £ + and £ - exhibit resonator-like behav- 
ior, i.e. they vary periodically in antiphase almost with- 
out damping with increasing interatomic separation. In 
Fig-ffl(c) is shown the dependence of £ ± (x = 0.32) [peak 
position of ^ in Fig. ^ (a)] on the interatomic separation 
for the atoms in the center of the metallic (9,0) CN. As 
the separation increases, £ + and £~ may differ greatly 
one from another in a periodic way. The maximal dif- 
ference is |f+ - 2(£- L + fH) « 4f- L - 10 8 , making 
the mixing coefficients C±(t) different and thus result- 
ing in a substantial degree of the entanglement of the two 
spatially separated atomic qubits. Figure [D (d) gives an 
example of the £ frequency behavior at x ~ x r = 0.32 
for the three interatomic distances [shown by the vertical 
lines in Fig. ^ (c)] corresponding to " antinode" relative 
positioning of the atoms in the center of the (9,0) CN. 

The ensuing spontaneous decay dynamics and atomic 
entanglement are presented in Fig.[21for the atoms in the 
center of the (9,0) CN. In (a),(b),(c) atom A is supposed 
to be initially excited while atom B is in its ground state. 
The entanglement is seen to reach the amount of 0.5 and 
to vary with time periodically without damping at least 
for the (reasonably long) times we restricted ourselves 
in our computations. As the interatomic separation in- 
creases, so the period of the entanglement oscillations 



does while no change occurs in the maximal entangle- 
ment. In (d) both of the atoms are initially maximally 
entangled [C± (0) = 1 while C T (0) = 0] and separated by 
the distance of 7.4 A. The entanglement is slightly larger 
in the case where C + (0) = 1, CL(0) = and no damping 
occurs as before. Note that the atoms can be separated 
by longer distances with roughly the same entanglement 
due to the £^ periodicity with interatomic distance. 

We have demonstrated sizable amounts of the two- 
qubit atomic entanglement in small-diameter metallic 
CNs. The entanglement greatly exceeds 0.35 that is 
known to be the maximal value one can achieve in the 
weak atom- field coupling regime [l3|. and persists with 
no damping for very long times. We attribute such a be- 
havior to the strong atom-field coupling and electronic 
structure resulting in the resonator-like distance depen- 
dence of the two-particle local DOS in metallic CNs. The 
scheme studied can be generalized to the multi-atom en- 
tanglement via the nearest neighbor pairwise quantum 
correlations, thus challenging novel applications of atom- 
ically doped CNs in quantum information science. 
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